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Abstract 

We discuss how in supersymmetric models with D and F-flat directions, a primordial monopole 
flux of order 10" 16 - 1(T 18 cm 2 sec 1 si 1 can coexist with the observed baryon asymmetry. A 
modified Affleck-Dine scenario yields the desired asymmetry if the monopoles are superheavy (~ 
10 13 — 10 18 GeV). For lighter monopoles with masses ~ 10 9 — 10 12 GeV, the baryon asymmetry 
can arise via TeV scale leptogenesis. 
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Unified theories of elementary particle interactions predict the existence of topo- 

logically stable magnetic monopoles with masses and magnetic charges that depend on 
the details of the underlying theory. In supersymmetric SU(5), for instance, the lightest 
monopole with mass of order Mq/c(g (— 5 x 10 17 GeV) carries one unit (27r/e) of Dirac 
magnetic charge (Here Mq ~ 2 x 10 16 GeV denotes the gauge coupling unification 
scale, and «g ~ 1/25 is the unified coupling constant at Mq). The monopole also has 
screened color agnatic fields Q.and the presence of baryon and lepton number violating 
superheavy gauge bosons leads to monopole catalysis of nucleon decay Q]. If the monopole 
arises from a partial unified model such as G422 = SU(4) C x SU(2)l x SU{2)r |l| (or 
G333 = SU(3) C x SU(3)l x SU(3)r), it carries two |(| (or three quanta of Dirac charge 1 . 
The absence in these models of gauge bosons which mediate proton decay means that, in 
principle, the symmetry breaking scale (M) can be considerably below Mq leading to 
lighter ( ~ 10 10 GeV) monopoles. In a supersymmetric framework intermediate scales are 
allowed provided proton decay via dimension five operators is adequately suppressed. Mod- 
els with intermediate scales have acquired renewed interest because they naturally ap pea r 
from compactification of superstring theories on Calabi-Yau manifolds and orbifolds 
from .D-brane constructions or from intersecting brane models jl^ . 

The fate of primordial monopoles is very closely linked to the history of the very early uni- 
verse. In hot big bang cosmology without inflation there is a serious cosmological monopole 
problem pointed out a long time ago ^jj. Here we list some proposed solutions: 

• An inflationary epoch [15I . ^| reduces the primordial monopole number density to 
utterly negligible values. The subsequent transition to radiation epoch may give rise 
to thermally produced monopoles, provided the reheat temperature is not too far 
below the monopole mass jl7|. 

• If the symmetry breaking which gives rise to monopoles experiences only the last 30 



or so e-foldings of inflation 







19j | . a measurable flux of monopoles comparable to or 



2l( may 



somewhat below the Parker bound [20j and the MACRO experiment limits 
)e present in our galaxy. A different scenario with similar objectives is discussed in 

3- 



1 Multiply charged monopoles also arise in superstring theories 
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Monopoles and antimonopoles get linked either by electromagnetic |23J or Z flux tubes 
24I ] and efficient annihilation occurs as a result of rapid contraction of these tubes. 
Another resolution may be non-restoration of grand unified symmetry at high temper- 
ature (25! ■ In [2^, it is proposed that if monopoles are produced in association with 
domain walls, the latter can sweep away monopoles and then d isap pear. A black hole 
solution of the cosmological monopole problem is discussed in \27\. 



Monopoles arise in supersymmetric models with D and F-flat directions 



which thermal inflation 
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31 
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33, 



34. 



28, 



111 



is followed by a huge release of entropy 



An initially large monopole number density can be reduced to the Parker 
bound or somewhat below it 



In this paper we wish to explore this last scenario, previously discussed in |29(, in which 
monopoles appear after thermal inflation is over and a huge amount of entropy is released. 
While this certainly helps with monopole dilution, our main challenge is to identify a frame- 
work in which an observable flux of primordial monopoles is compatible with the observed 
baryon asymmetry. This is particularly challenging for superheavy monopoles with a cor- 
respondingly large symmetry breaking scale, so that the final temperature after thermal 
inflation can be quite low, of order a few MeV. Following 37], we employ a modified Affleck- 
Dine (AD) scenario to realize the observed baryon asymmetry. Although our discussion is 
mainly focused on monopoles with mass ~ 10 14 — 10 15 GeV, it can be adapted, as we will 
show, to monopoles that are somewhat heavier (~ 5 x 10 17 GeV) or significantly lighter 
(~ 10 9 - 10 13 GeV) Q- 

Consider the superpotential 



W inf D A 



Mf 



2n-3 



+ /3H u H d 



M, 



2m-l ' 



(1) 



where the scalar components of and 0, called flatons [3jj, acquire non-zero vevs (M) and 
break the underlying gauge symmetry G. M* is the cutoff scale, H u d are the MSSM Higgs 
doublets, and n, m are integers suitably chosen to yield the desired symmetry breaking scale 



M and the MSSM [i term respectively 



28, 



37 



39 
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4l| . The second term in Eq. (P) also 



plays an important role in reheating after thermal inflation. 

The zero temperature effective potential of flatons (we use the same notation for the 
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superfield and its scalar component) is given by 



V{<f>) = £ - M] 



+ 



A0 n 



■n — 1 



Ml 



(2n-3) 



+ 7T 



A0™~ i 



Mi 



(2n-3) 



Mi 



2n-3 



+ c.c, (2) 



where \i\ is the false vacuum energy density, such that V(M) = at |(0)| = M, M s is 
the soft supersymmetry breaking mass parameter, and \A\\ < M s . Minimizing the effective 
potential along the D-Hat direction |(0)| = (0)1 yields the symmetry breaking scale 2 , 

N 1 l/2(n-l) 

\A X \ + ^\A X \ 2 + 4(2n-l)M s 2 

l/2(n-l) 

(3) 



M 



(2n-3) 



M, 

2(2n - l)n\ 



M.Mi 



(2n-3) 



v/(2r2 - l)nA 

where n > 2 and for simplicity, we assume \A\\ < M s . Some typical values of M s , M* and 
M that we will consider in this paper are listed in Table I. For non-zero temperature T the 
effective potential gets an additional contribution given by 



v T {4>) 



(^) E(- 1 ) F /°° rfxx2ln ( x - (-!) Fex p{- 



v 2 , MM 

X + ^2 



T 1/2 



} 



(4) 



where the sum is over all helicity states, (— 1) F is ±1 for bosonic and fermionic states, respec- 
tively, and Mj(0) is the field-dependent mass of the zth state. For T the temperature- 
dependent mass term is a T 2 |0| , where a ~ 0.1. For T > T c = M s /y/a the potential 



^(0) = + 2(-M 2 + ( xT 2 )|0| 2 + 2n 2 A 2 ^-^ 



+ ^aA- 



r2n-3 



+ C.C 



(5) 



develops a minimum at |0| = 0, the gauge group is unbroken and hence there are no 
monopoles. The false vacuum energy density /i 4 , (~ M 2 M 2 ) drives thermal inflation and the 
universe experiences roughly ln(/i /T c ) e- foldings (~ 12 for M ~ 10 14 GeV). 

As the temperature T falls below the critical value T ~ T c , the mass-squared term for 
turns negative, and rolls from the origin to M, thereby ending thermal inflation. Monopoles 
are expected to arise as a consequence of symmetry breaking through the Kibble Mechanism 
The field performs damped oscillations about the minimum at M and subsequently 
decays. During these oscillations, the universe is matter dominated with energy density 



Here the phases e, a and a of A\, <f> and <fi are taken to satisfy the relation e + na + na = it. 



4 



p ~ M s 2 (0) 2 (t 2 /t 2 ), where t c represents the cosmic time at the phase transition. A large 
amount of entropy is released by the decay of field which helps dilute the initial monopole 
density |29j. To estimate the initial monopole number density um one usually makes the 
plausible assumption that a correlation size volume ~ £ 3 (~ T~ 3 ) contains on the order of 
one monopole QQ Th en Q,Q 

P 



n M 



10- 2 , (6) 



LT 3 J initial (4tt/3)£ 3 T 3 

where P is a geometric factor of order 1/10. It has been suggested |44j that £ may be 
somewhat larger than T c _1 , in which case fewer monopoles would be produced. A more 
drastic reduction in the initial monopole number density is achieved by assuming that a 
single monopole is produced per horizon volume 4^|- Since the horizon size during monopole 
production is larger than T~ x by a factor Mp/M, a large suppression (by a factor of order 
(M/Mp) 3 ) of the monopole number density becomes possible. We will not use this last 
suppression mechanism for superheavy monopoles, but will keep in mind the fact that the 
initial monopole number density estimates can have large uncertainties in them. 

One may ask whether r in can be reduced through monopole-antimonopole annihilation 
14,4^. In the context of flaton models this has been investigated in j3] which reached the 



conclusion that there is no significant annihilation during the epoch of domination. 

We now consider monopole dilution through entropy production from the decay of the 
flaton field. If the entropy increases by a factor A, the final number density of monopoles is 
given by 

= P M 1 - r \-l 9*i T final) m 

fmal ~ IT 3 \ final in g*{T c ) ' [) 
where g*(T) represents the degrees of freedom at T. The flaton decay proceeds predom- 
inantly via the coupling (3H u H d ((j)(f)) m / ' M^ m ~ l . The decay width ((f) — > Hj) is given 
by 

1 f M\ A{2m - l) M 2 



where = 2\j2{n — 1) M s is the flaton mass. The final temperature after thermal inflation 
can be expressed as 

. - „, (M X 2 ' 2 "- 1 ) M /Mp 

zv* 0.3 0.036/J" ( W J ^TT^V^' (9) 
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where Mp = 2.4 x 10 18 GeV is the reduced Planck mass. Note that the last term in Eq. 
generates the effective MSSM \i term f3{M/M*) 2m - l M. Eq. © can be rewritten as 

/ u 2 \ 1 IWf 

T f ~ 0.036 ' 



M) ( n -if 4 V M s ' 

1/2 



/ n \2 /10 14 GeV\ /950GeV\ v ' Ar , p . , 

aM (iiif) HH HH ^ (f«n = 4;m = 3). (10) 

The entropy released by ffaton decay is estimated to be 

A ~ M 

9*{T c )T*T f ' 

where we have used Eq. (fTUj) and g*(T c ) ~ 200. Using Eqs. (|TT|) and (0) we get 
« m-28 f r - A /0.1\ 3/2 /KPGeV\ 3 / /i \* / \ 1/2 

r ^" 6x10 (ifOItJ (, — ivZ — J U (gsoGevJ ' (12) 



for <fc(7»~10 |47|. 

From Eq. (jlOj) we note that if M increases above 10 15 GeV, the \x parameter and M s also 
have to increase (see Fig. in order that the final temperature Tf stays above 10 MeV for 
successful nucleosynthesis. Thus, some fine tuning of the soft supersymmetry breaking higgs 
parameters will be required to implement electroweak breaking. Monopoles with masses 
greater than or of order 10 17 GeV can gravitationally clump and will be considered later. 

We should make sure that the final monopole number density is consistent with the 



MACRO bound [2]J which is more stringent than the well known Parker bound |48|, |4£ 
For superheavy magnetic monopoles (10 11 GeV < tjim ^$ 10 17 GeV) moving with speed 
v M ~ 3 x 10- 3 c(10 16 GeV/m M ) 1/2 HQ, 

the MACRO limit on the flux Fm is 

F M < lO-^cm^sec^sr- 1 , (13) 

which corresponds to r fina i < 2 x 10~ 26 (3 x 10~ 3 c/v M )(m M /10 16 GeV) 1 / 2 . For n = 4, m = 3 
in Eq. ([!} and with the remaining parameters given in Table I (corresponding to M — 10 14 
GeV), an entropy release factor A of order 10 22 saturates the bound in Eq. ([13)1. with 
r in ~ 10~ 2 (see Eq. The variation of M s and M with the initial monopole number 

density given by Eq. (jSJ, such that the monopole flux bound (Fm — 10 _17 cm _2 sec~ 1 sr~ 1 ) 
is saturated, is shown in Fig. [2(a). For Fm = 10 _18 cm _2 sec _1 sr _1 M s vs. M is plotted in 
Fig. Mb)- 
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10 14 


950 


10 16 


0.01 


10 3 


0.01 


0.1 


lO" 2 


lO" 17 


10 15 


3000 


6.4 x 10 16 


0.3 


1.8 x 10 4 


0.0003 


0.02 


10"2 


1Q -18 


10 16 


3000 


10 18 


0.03 


2 x 10 4 


0.0003 


0.02 


10~ 2 


1Q -16 


2 x 10 16 


2500 


2.4 x 10 18 


0.01 


1.6 x 10 4 


0.0003 


0.02 


10~ 2 


io- 17 



TABLE I: A set of parameter values for which the expected monopole flux is at or below the 
MACRO bound. 

The release of such a large amount of entropy (Eq. (fTTj) ) certainly washes away any 
pre-existing baryon asymmetry. Also, for M > 10 13 — 10 14 GeV the final temperature 
is quite low 3 {Tf ~ 0.01 GeV), so that the sphalerons are ineffective, and the standard 
leptogenesis scenario [3] does not apply. A different mechanism for generating the desired 
baryon asymmetry must then be found. This problem has arisen before and discussed by 

for instance, new particles beyond the MSSM spectrum 



several authors 



33 



53 



are considered whose out of equilibrium decay can give rise to the baryon asymmetry. Here 
we rely on a modification of the AD scenario proposed in J^tJ (and subsequently in 5(]]) 
in which a dilution factor A -1 ~ 10~ 17 — 10~ 18 and symmetry breaking scale M ~ 10 10 ~ n 
GeV are considered 4 . In our case, as we saw in Eq. (jllj) . the monopole problem requires an 
even greater amount of entropy production, especially if ri n ~ 10~ 2 . 

To implement the scenario discussed in |37J in our case, we couple <p to the right-handed 
neutrino superfields N{. In the MSSM notation, consider the couplings 5 



W 



W inf + YiLH d e + Y D LH U N + X^N 2 + 



(14) 



where we have suppressed all generation and group indices, and the ellipsis represent terms 
in the superpotential which will not participate in the analysis. We assume that both during 
and after thermal inflation, the squark fields do not acquire non-zero vevs. 



3 Tf as low as this can lead to conflict with LSP cosmology [5(J • One way out may be to introduce an axino 
LSP 0. 

4 For noncquilibrium effects in AD mechanism, see |57| . 

5 In practice, the superpotential is invariant under G, and Ni belong in some appropriate G representation. 
Later we will provide an example for G and discuss how the terms in Eq. 114(1 arise. 
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Consider the zero temperature F-term potential 



V F = \Y D H U N + Y { E d ef + \Y D LH U + 2A^iV| 2 + 



YiLe + (3H U 



(00) : 



+ 



+ 



-r\3 



Y D LN + (3H d - 



MP 



4A 



-A 3 



"T\2 



Mi 5 



3f3H u H d 
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3 / ,-rv2 
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+ A A^ + 3/3H u H d 



(15) 



* * 

where we have used the same notation for superfields and their scalar components and set 
n = 4,m = 3 in Eq. (Q). The D-terms (Vd) and the supersymmetry breaking parts of the 
potential {%^ S y) are 

V D + V sisy = g 2 (\H U \ 2 - \H d \ 2 - \L\ 2 ) + m\ |L| 2 - m\ u \H U \ 2 - 2M 2 |0| 2 



+ A D Y D LH U N + A A . A^0iV 2 + A (3H u H d 



-r\ 4 



MS 



+ ^aA 



Mi 5 



+ ••• +c.c, (16) 



where the soft masses and A-terms are all taken to be of order M s , and D-flatness along 0, </> 
direction has been taken into account. For simplicity we will confine ourselves to a single 
generation picture. 

We assume that initially all fields are held at zero during thermal inflation. To implement 
the scenario, the AD-field, parameterized along the D-flat direction by LH U = ip 2 /2, acquires 
a vev ~ m\j \ Yd\ 2 provided we set m 2 ^ = {m 2 Hu — m 2 L )/2 > (0 is still at zero). Assuming 
this happens before the flaton field starts to roll down, the term Ae>YdLH u N + c.c. induces 
a vev for N which, in turn, triggers the slow roll of towards its minimum. As <ft increases, 
N follows its instantaneous minimum 

Y D LH U 



N 



2X (j> (p 



(17) 



Using Eq. (fTTj) the effective potential reduces to 6 



V = V F + V D + V sisy 

2 



4A 



M5 



1 14 



2Mt 







ill 5 



1 12 



+ 



4A 



+ 



A, 



-A D + 2X* 



3 -T* 



M*</> / 8A^ 



+ ^4aA 



M 5 



+ c.c 



(18) 



For simplicity e and Hd are assumed to be at their zero location. 
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It follows, after some algebra, that for <C M, ip is given by 



and its phase is related with that of 4> from the term 

^Sf)+- (») 

Note that at this stage the third term in Eq. (|20|) is negligible. 

When (J) becomes comparable to M (with the temperature still around T c ), it yields a 
positive mass squared term for if), provided that 

\f>%-<>»- < 21 > 

This helps the LH U direction to reach its true minimum at zero 7 . Simultaneously, the phase 
of "0 gets an important contribution from the third term inside the bracket of Eq. ([20)1 . This 
change then initiates an angular momentum in LH U which turns out to be the amount of 
lepton asymmetry (ri£ = i(ifj*ip — ipip*)) generated from the AD mechanism. The difference 
between (A\./2M a — Ao/M s ) and 2A* would be the source of CP- violation. 

To protect til from being erased by potential lepton-number violating processes, in Ref. 
] it is assumed that the AD field should completely decay while the universe is still 
dominated by the energy of the flaton field. We assume that the decay width of the AD 
field is such that the electroweak sphalerons 58] can convert a fraction of til into tib, the 
baryon number density. 

The decay of <fi yields a final temperature of around 0.01 GeV. The release of entropy (es- 
timated earlier) responsible for the dilution of monopoles also dilutes the baryon asymmetry. 
The final baryon asymmetry is approximately given by 



S 3 Si n S final 

~^ Tj ~tAb*Tt?WL t (22) 



7 From the tabulated values of A, f3 (corresponding to M = 10 14 GeV) and Eq. (1221) we notice that this 
happens before <f> actually reaches its minimum (as A <C (3). 
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where we have used Eq. f)19|) . The seesaw relation for the light neutrino mass is 8 

y2 2 

m * - tu- (23) 

where we employ a basis in which the light neutrino mass matrix is diagonal and, for sim- 
plicity, only a single LH U flat direction is assumed. Substituting Eq. in Eq. 

s 3 V 3 m v M 2 M s ' 

where n L /n^, ~ 0(1) is assumed. 

To obtain the required Ub/s, the neutrino mass m u in Eq. (J24j) turns out to be several 
orders of magnitude below the scale for atmospheric and solar neutrino oscillations. From 
Eq. 



a b i /! v 2 in,.'!) 



n B 
s 



iU ll74GevJ V mu )\ M ) ^O.OlGeVy'' [Zb) 



7 — 



where we have set ~ M s , and from the observed baryon to photon ratio ri^/n, 
(6.0965 ± 0.2055) x 10" 10 H, n B /s ~ (n B /n 7 )/7.04 In Fig. Hwe display the allowed 
region for m v , with the range of M and Tf restricted by the monopole flux corresponding 
to Fig. |2](a). 

In order to make the discussion more explicit and to realize a scenario with superheavy 
symmetry breaking scale corresponding to n = 4, m = 3, we consider a realistic model with 
gauge symmetry G422 = SU(A) C x SU(2) L x SU(2)r |l|, whose breaking yields monopoles 
which carry two quanta of Dirac magnetic charge [6j]. By the same token, in principle, there 
also can exist color sing let states ^ ch arg e ± 7/ 2 fl We have 

as low as 10 14 GeV is compatible with proton lifetime limits. Indeed, for M ~ 10 14 GeV, 
we estimate the proton lifetime to be of order 10 34 yrs, taking into account the operators 
discussed in |8JJ. The quarks and leptons are unified in the representations Fj = (4,2, 1)^; 
Fi = (4, 1, 2)j of G422, where the subscript i (= 1, 2, 3) denotes the family index. The Higgs 
sector consists of 

h = (1,2,2); 0= (4,1,2); 0= (4,1,2). (26) 



From Eq. 114(1 . it is seen that the effective superpotential responsible for generating the neutrino mass is 
W u = -(Y D LH U ) 2 /(2A 0) after integrating out N. 

If G422 is replaced by G333 = SU(3) C x SU(3)l x SU(3)r we find monopoles carrying three quanta of 
Dirac charge 
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TABLE II: Discrete charges of various superfields. 



The superpotential W is given by 
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-FF 
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+ «P» + <>*>** + jr™ [ft 



(27) 



where a discrete symmetry Z4 x (see Table II) has been introduced to realize the desired 
D and F-flat direction 10 . The color sextet superfield T> = (6, 1, 1) = (D C ,D C ), where D c = 
(3, 1, 1/3) and D° = (3, 1, —1/3), is introduced to provide heavy mass to the components 
d c H ,d c H of 0,0 0,Q- 

The interaction of with iV in this particular example is determined by the term (70 2 + 
K(f) 2 )N 2 /M# in the superpotential instead of X^N 2 in Eq. (|14j). Thus as long as <fi is held at 
zero, N does not acquire any vev through the term AYFhF and from Vp even when the AD 
field has a vev 11 . The field starts to roll down when T ~ M s as discussed in the context of 
Eq. jnj). As increases, would receive a vev ~ (y(-0) 2 /4(7 + k) |0|)(M*/ |0|) similar to 
Eq. (HZD, with A ~ (7 + k)(|0| /M,) = 7 '(|0| /M„). The relation in Eq. (H| now becomes 



"kin 



M S M 
M ± 



101 • 



(28) 



The rest of the discussion for an estimate of the observed baryon asymmetry is very similar 
to what we already have before Eq. (|25j) . 



The baryogenesis scenario we have discussed, following 37[, can also be employed for 
monopoles with masses ~ 5 x 10 17 GeV, corresponding to GUT symmetry breaking scale 
M ~ 2 x 10 16 GeV. The gravitational force on such monopoles exceeds the magnetic force 
and it seems plausible that they would clump in the galaxy, and perhaps even contribute 



10 To avoid topologically stable domain walls we will assume that the discrete symmetry is explicitly broken 
by higher dimensional operators. 
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The AD field, tp, in this example is a flat direction chosen along the neutral components of F and h (H°). 
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to the dark matter in the universe. (This may not be plausible for monopoles that catalyze 
nucleon decay). However, this latter possibility is disfavored by the MACRO bound [9]. To 
see this, assume that the monopole energy density pu < Pb, the baryon energy density. 
Taking a local density enhancement factor of order 10 5 , the local flux is estimated to be 
(with v M = 10~ 3 c) 

F M < 6 x 1(T 13 ( 1Ql7GeV ^ cm^sec^sr- 1 ; m M ~ 10 17 - 10 18 GeV , (29) 

V m M J 

which is in strong disagreement with the MACRO bound. Thus, pu 10" 4 Pb, to be 
consistent with the MACRO bound. To achieve a monopole flux < 10 _17 cm _2 sec _1 sr _1 , we 
estimate that rji na i < 10 -31 . The parameters chosen to achieve this are displayed in Table 
I. 

Let us now consider 'lighter' monopoles with masses of order 10 9 — 10 10 GeV, correspond- 
ing to symmetry breaking scales ~ 10 8 — 10 9 GeV. (In models with D and F-flat directions 
the symmetry breaking scale is expected to be > 10 8 GeV). In this case the amount of 
entropy released following thermal inflation is considerably smaller. Namely, from Eq. (fTTj) . 
A ~ 10 6 , so that a monopole flux of order 10 -17 cm _2 sec~ 1 sr~ 1 (this is consistent with 
present bounds from MACRO 



2l|, SLIM |63( and AMANDA |64| experiments) requires 
that the initial relative monopole number density r in must be sufficiently small, as shown 
in Fig. E] (see discussion following Eq. (jUJ)). Following j^J, the observed baryon asym- 
metry can now be explained via TeV scale leptogenesis, since the final temperature Tf ~ 
TeV. An initially large lepton asymmetry can survive the moderate amount of dilution so 
as to produce a final baryon asymmetry consistent with observations. Finally, we see from 
Fig. |H that a flux bound of intermediate mass monopoles (~ 10 11 — 10 13 GeV) requires 
that r in ~ 10~ 17 — 10 -7 . The baryon asymmetry here can be achieved either through TeV 
scale leptogenesis j^] or the original AD scenario [3| where Tf is of order 100 GeV. For 
monopoles of mass ~ 10 13 — 10 14 GeV, the baryon asymmetry can be generated through the 
modified Affleck-Dine mechanism. 

To summarize, magnetic monopoles appear in a variety of unified gauge models with 
a wide range of masses. In models of thermal inflation, monopoles with masses of order 
10 9 — 10 18 GeV appear when the cosmic temperature is in the TeV range. Their subsequent 
dilution through the release of huge amount of entropy poses a challenge for baryogenesis. 
We have discussed a class of realistic models in which the observed baryon asymmetry can 
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co-exist with a primordial monopole flux which can be detected with large scale detectors 
such as ICE CUBE |s5 1 . 
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FIG. 1: M s vs M (for superheavy monopoles), with final monopole flux Fm = 10 cm 2 sec x sr . 
Here 100 GeV < p < 20 TeV. 
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(a) u. < 2 TeV 



(b) n < 2 TeV 
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FIG. 2: M s versus M with A = 1CT 2 and 10~ 3 < f3 < 1. Here r in = 10~ 2 . The monopole flux 
Fm = 10 _17 cm _2 sec~ 1 sr _1 for plot (a) and 10~ 18 cm _2 sec _1 sr _1 for plot (b). 
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FIG. 3: m v versus M for Fm = 10 _17 cm _2 sec~ 1 sr~ 1 and rj„ = 1CP 2 (corresponding to Fig. |U(a)), 
with n B /s ~ (0.837 - 0.895) x lO" 10 59]. 
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FIG. 4: Ti n vs. M (for relativistic monopoles), with Fm — 10 17 cm 2 sec ^sr \ 100 GeV < M s < 
1 TeV and 100 GeV < fi < 2 TeV. ru represents the relative monopole number density assuming 
production of one monopole per horizon volume. 
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